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t-H Abstract. The Hirota equation is better than the nonlinear Schrodinger equation when ap- 

proximating deep ocean waves. In this paper, high-order rational solutions for the Hirota 
equation are constructed based on the parameterized Darboux transformation. Several types 
of this kind of solutions are classified by their structures. 
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1. Introduction 

In recent decades, the rogue waves fl] (or freak waves (21), have gained more and more interests 
in the study of physics. One of the reason is that such kind of waves emerges spontaneously 
and frequently in the ocean without any sign, "appears from nowhere and disappears without 
a trace J3]". Recently, besides oceanography, the rogue waves are also found in the optical 
systems when generating highly energetic optical pulses |4j[5], as well as in the superfluids and 
femtosecond pulse (6 10 , etc. The growing consensus, which is from the oceanic and optical 



\Q rogue waves, arises from to the modulation instability (MI) [II 15 which is a key property 
for many integrable equations. One of the typical example is the nonlinear Schrodinger(NLS) 
equation 

O 

m i qt + 2\q\ 2 q + q xx = 0. (1.1) 



For this equation, the lower-order rogue wave has been obtained 16 with a good experimen- 
tal simulation. People are more interested in the higher order waves, because its numerous 
and complicated properties. Some special higher order rogue waves have been constructed by 



c3 Darboux transformation 17,18 , and some more general cases (multi- Peregrine) were well stud- 



ied 19-23 . Unlike soliton solutions siting on the zero background, these higher order rogue 



waves not only sit on the non-zero constant background, but lots of behaviors are quite different 



at points on the space-time plane 23 . On the other hand, the rogue waves are the limit of 



a sequence of the breather solutions 18 , 24 . These new studies about rogue waves enrich the 

conception of rogue waves and lead to understand further towards this mysterious phenomenon. 

However, the phenomena of the ocean are significantly complicated, such as the depth of the 

sea, bottom friction, viscosity, etc, these factors all need to be taken into account when we do 



the simulation. Therefore we need to add higher order terms and other nonlinear effects 25,26 
The Hirota equation can be viewed as a generalization for the NLS equation, when the higher 
dispersion and time-delay changes is taken into account. 
In dimensionless form, the Hirota equation is given by 

iqt + a(2\q\ 2 q + q xx ) + i(3(q xxx + 6\q\ 2 q x ) = 0, a, (3 G E. (1.2) 
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The variable t is the propagation variable and x is the retarded time variable in a moving frame. 



The phase function q is the envelope of the wave field. The two terms in (1.2) that enter with a 
real coefficient j3 represent the third-order dispersion (q xxx ) and a time-delay correction to the 
cubic term (Iqi 2 ^), respectively. We can see that 



2, 



(1) a = 1,(3 = 0, (1.2) becomes the NLS equation: iq t + 2\q\ 2 q + q xx = 

(2) a = 0, /3 = 1, (1.2) becomes the modified KdV (mKdV) equation: 
Qt + Qxxx + 6|g 



% = 0. 

In some cases, if an integrable equation is a combination of two known equations that are 
also integrable, the low-order solutions can be derived from solutions of those two equations. 
But for high-order rogue waves, it does not work. Therefore, we provide the parameterized 
Darboux transformation to derive the high-order rational solutions. 

The organization of this paper is as follows. In section 2, based on the parameterized Darboux 
transformation, the general formation of the solution is given in terms of the Gram determi- 
nants. In section 3, we construct the higher order rogue waves from a periodic seed with a 
constant amplitude and make a classification based on their structures. We provide further 
discussion in section 4. 

2. The Darboux Transformation for the Hirota equation 
In this section, we give the explicit Darboux transformation for the Hirota equation. 



The Lax pair for the Hirota equation is as follows 27 



Mip 
Nil), 



(2.1) 



where 



'I \ U ._. I ~^ <? \ A' _ C V! f ' \2 

r iX 



V=( £ ),M=( J* ,\ ),N = V 3 X i + V 2 \ 2 + V 1 \ + V , (2.2) 



/ -Ai(3 \ / -2ia Aq(3 \ _ ( -2iqr(3 2iq x (3 + 2qa \ . , 

3 V Aif5 ) ' 2 V A ^ 2ia J \ - 2ir *P + 2ra 2i( F$ ) ' 

y _( -iqra + (3(-qr x + q x r) iq x a - (3(q xx - 2q 2 r) \ _ ^ 

y -ir x a + /3(-r xx + 2qr 2 ) iqra - /3(-qr x + q x r) ) ' 

Here A is a spectral parameter with complex value, and ip is the eigenfunction corresponding 
to the A of the Hirota equation. 

In terms of zero-curve equation M t — N x + [M, N] = 0, equations for q and r are as follows: 



q t = ia(2\q\ 2 q + q xx ) - (3(q xxx + 6|g| 2 g x ), 
r t = -ia(2\r\ 2 r + r xx ) - (3(r xxx + Q\r\ 2 r x ), 



These two equations admit the constraint 

r = -q*, (2.6) 



here * denotes the complex conjugate. After the constraint, the system (2.5) becomes one 
equation, which is exactly the Hirota equation. 
Next, we will construct a gauge transformation 

^ [1] = TV, (2-7) 
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such that, after this gauge transformation, 



-- (T x + TM)T-\ 
(T t + TN)T-\ 



(2.8) 



is solvable for tp'- 1 ' and Af' 1 ', N^> are of the same form as M, N. Moreover, we want M M, N^ 1 ' 
also admits the constraint (2.6). Therefore, we can get new solution of the Hirota equation 
from one solution. 

Let {ipj',1 < j < n} be eigenfunctions for the Lax pair (2.1) which are corresponding to 
eigenvalues Xj,_ 



>J = 1,2, 



, n: 



^3 



<j>j2{x,t, Xj) 



J = 1,2, 



.71. 



(2.9) 



Moreover, if ip 3 - is a eigenfunction for eigenvalue Xj, then 



#:= 



«5i 



is the eigenfunction corresponding to A*. 

Choose a sequence of eigenvalues { Xj\j = 1, 2, 3, ..., 2n}, such that A^-i = A 2n and {tpi,ip2i "03, 
is the corresponding eigenfunctions. We use a modified mechanism described in [18] to construct 
the parameterized Darboux transformation for the Hirota equation: 
Theorem 1. Given Ai G C, and 



01 



h2 



is an eigenfunction corresponding to Xx of (2.1), denote 
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with A 2 



•'ii 
hi 



h-2 



then 



'0 [1] = T^i 



is again a solution of (2.1) and 



qi = q + 2% 
is a new solution for the Hirota equation. 



A, 



Ai0n 

A2021 



i>2i 



(2.11) 



(2.12) 



This formula is of the same form as in 18 . For high-order rational solutions q n ,r n , calculation 
is quite cumbersome, but one can still check these equations with the aid of modern tools such 
as "maple" or equivalent, and also numerically. 

Theorem 2 (N-fold). By applying the gauge transformation repeatedly for N times, choose 
eigenvalues {Ai, A 2 , A 3 , ..., A 2n }, such that A^i = A 2n , n = 1,2, ...,N, and ipj is an eigenfunc- 
tion for Xj. Then we can get N-fold solution for the Hirota equation: 



Qn = q 
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(2.13) 
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As a consequence, high-order rogue waves for the NLS equation and the mKdV equation can 
be derived. 

3. The High-order rogue waves 

In this section, we give high-order rogue waves for the Hirota equation from non-zero seeds. 
By using the principle of the superposition of the linear differential equation, we can get the 
new eigenfunction corresponding to Aj: 



i>i 



<Pji 

Vj2 



fji = dt4>ji + d 2 <pj2, <Pj2 = -di4>* 2 + d 2 (j) 



with 
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(3.1) 



a 2 + 4c 2 + 4A,-a + 4A 2 , c 2 



a\j + 2/3A 2 — -aa — 8 c 2 + -8 a 2 — XmB. 



Here a, 6, c, s,- 



j. -- -v ■ -r---j 2 
0, 1, 2, . . . are the constants, e is infinitesimal, -01 



of (2.1) for the non-zero seed q = ce tp } p = ax + bt, b = a(2c 2 



[0ii ; 0i2) T is the solution 
a 2 ) + /3(a 3 -6ac 2 ). 
It is well known that high-order rogue waves can be obtained by the Taylor expansion of the 
breather solutions. In [18], the scheme for nonlinear schrodinger equation is established. Due 
to the complication and tediousness of the explicit formula, we describe the high-order rogue 
waves by their structures, explicit formula for solution qi and 52 will be given in the Appendix. 
Moreover, in the rest of this paper, we will dedicate to the classification of rational solutions. 
In our classification, we can see that by choosing different values for parameters d\ and d>2 
different types of solutions are obtained. 

(1) When d\ = d 2 = 1 (after expansion), rational solutions of any order n have similar 
structure. 

(a) There are ^^ — - — 1 local maximum in each side of the line t = 0. 

(b) Starting from —00, before the central high amplitude (the global maximum), there 
is a sequence of peaks with increasing height. In detail, the number of first peaks 
is n, then there is a row of n — 1 peaks and so on. 

4 



(c) The structure is symmetric with respect to time t. 



We called this type of rogue wave the fundamental pattern (same as in 18 Figure 1) 
(2) Ring structure. 



In (3.1 ), when n > 3 and the principle coefficient for order n rational solution s n _i 3> 
1, while other coefficients Si are all zero. The structure for rational solutions consists of 
two parts. The outer circular shell consists of 2n — 1 first order rational solutions, while 
the center is an order (n — 2) rational solution of fundamental pattern. This structure 



is similar to the "wave clusters" in 28 . Figure 2 give a visual description for the ring 
structure. 
(3) Triangular structure. 

Another typical type of rational solutions can be obtained by choosing the first non- 
trivial coefficient si ^> 1, while others are all zero. The order n rational solutions have 
a structure of equilateral triangular type with n peaks in each edge. Each single peak 
is of the same height and the total number of peaks is n ' n 2 + ' (Figure 3). 

Remark 3.1. When n = 2, the ring structure and triangular structure have the same shape: 
there is a single triangular with 3 = 3n — 3 peaks in the triangular structure, and in the ring 
structure, since n — 2 = 0, there is no basic mode in the middle and the outer ring consists of 
3 = 2n — 1 peaks. 

These three types of structures listed above are the basic scheme for high-order rogue waves. 
Another interesting fact worth emphasized here is that: high-order rogue waves are usually the 
combination of these three types: 'fundamental pattern", "ring" structure and ' 'triangular" 
structure. In the rest part of this paper, we will discuss these combination types and their 
dependence on the free parameters. 

(1) As we mentioned before, the principle coefficient "s n _i" determined the ring structure. 
Based on this, if only the preceeding coefficient "s n _2" is non zero, the high-order rogue 
waves have a double-ring structure. More explicitly, there are two rings with 2n — 3 
single peaks in each and one (n — 4) basic mode in the middle (Figure 4). For example, 
(c) shows this type of combination for \q§ | 2 . There are two outer ring structures with one 
2-nd order basic mode in the center (Figure 1 (b)), and each individual ring is consisted 
of 9 single peaks. When n — 4, there is no fundamental pattern since n — 4 = (a). 

(2) Triple-ring structure is determined by the coefficient s n _3. When s n _3 ^ 0, the order n 
rogue wave has a triple ring structure with n single peaks in the outer circular and an 
order n — 1 double ring structure in the middle (Figure 5). 

(3) If S\ and s n _i 7^ 0, the phase space of high-order waves can be decomposed into two 
parts. The outer circular shell has 2n — l single peaks while the center is an order n — 2 
rogue wave with triangular structure (Figure 6). 

(4) If both s n _i and s n _2 are non zero, then the order n rogue wave is a composition of 
order n and order n — 2 ring structure such that the outer cell has 2n — 1 single peaks. 
And the inside is an order (n — 2) ring structures (Figure 7). Here we need to point out 
that although the shape is similar to double-ring structure described above, they are 
of different pattern in combination. In double-ring structure, both two outer circulars 
have 2n — 3 single peaks in each. While here, the outer cells have 2n — 1 and 2n — 5 
single peaks in each respectively. 

(5) Other parameters also play an important role in determining the structure of high- 
order rogue waves. In other words, by choosing the parameter properly, we can get 
more complicated structures. Here we give one more example. It is a composition of 

5 



order n ring structure and n — 2 double-ring structure (Figure 8). Moreover, it is a 
complete decomposition in the ring structure. That is, each local maximum is of the 
same value and there is no fundamental pattern. 

Remark 3.2. There are several coincides between different patterns, that is due to the special 
value of n. For example, when n = 4, then n — 3 = 1, therefore, the triple-ring structure of the 
same shape as ring structure is for triangular structure. 

Beside different structures, there are some properties that are preserved by all types of rogue 
wave solutions. First, the total number of peaks in order n solutions is ^^ — -. Secondly, the 
structures depend on every parameter. The principle coefficient s n _i leads to the ring structure. 
The first nontrivial coefficient Si leads to triangular solutions. They are very similar to the 
rational solutions of the NLS equation apart from the finite tilt with regard to the axes. We 
find that the main effect of the third-order dispersion and a time-delay correction terms added 



to the NLS equation is to introduce a "tilt" to the NLS equation solutions 29 , namely, the 
tilt is related to the finite value of parameter /3, because if we set /3 = the tilt is zero. We 
can shorten the length of the "ridge" by taking the value of a properly, moreover, of interest 
it is to find that on the central line where t = the shape is independent of /3 and a. That 
is, the adding the mKdV terms to the NLS equation do not have any impact in the figures of 
the rogue waves at the t = where the figures have its maximum. Consequently, the overall 
maximum value of the amplitude figure is qj(0, 0) for any value of j3, which is indeed the same 
as the NLS equation. 



4. Conclusions and discussions 

In this paper, we use the parameterized Darboux transformation to construct the high-order 
rational solutions for the Hirota equation. Based on the theoretical result, we classify the high- 
order rogue waves with respect to their intrinsic structure. We find out that there are three 
principle types, basic mode, ring structure, and triangular structure. Moreover, the composition 
of the three principle structures are held by high-order rogue waves. As we mentioned in the very 
beginning, the Hirota equation is a generalized NLS equation. The ring structures obtained in 



this paper are similar to the "atom" in Ref 28 . This explains the generalization in terms of the 
solution. On the other hand, by changing the free parameters in the Darboux transformation 
we can get more complicated structures. It would be promising to ask whether there are some 
analogues between high-order rogue waves of different integrable equations? Also, it is essential 
to find more conserved properties for such kind of solutions. All of these will help us in a better 
understanding in the deep ocean wave. 

On the other hand, the most remarkable feature of this equation is that it includes both the 
NLS equation by restricting the parameter = and the mKdV equation by restricting the 
parameter a = 0. Namely, by taking the different values of the parameters, we can also obtain 
high-order rogue waves of the NLS equation and mKdV equation, with regard to this point, we 
shall give the detailed describing in other papers respectively. 
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5. Appendix 

The formulas for g 1; q 2 are given under the condition Sq = 0,a = 10, P = l,a = 0,c = 1 and 
a — 5, j3 — 1, «o = 0, a — 3.88, c = 1 respectively. 

an,-, 3 - 1744 t 2 + 48 tx + 160 it -Ax 2 
qi = e MU - 



l + 4x 2 + 1744£ 2 -48£x 



i(3.88 rr-30. 140928 t)^l 

Ai =45 + 351839417.0 £ 6 - 252396.4448 t 2 x 2 + 33987.4508 £ 2 x 4 
+ 139.4688 x 5 £ + 4344315.56 t 5 x + 49254.311 x 3 £ 3 
+ 1326.7968 xH + 5998449.26 tV - 996552.958 xt 3 
+ 3325.2480 tx + 1307537.37 it A x - 127944.917 it 2 x 

- 4780.80 it + 7408.583 it 2 x 3 + 1802722.57 it 3 x 2 
+ 144 i Sl - 180 x 2 - 72088.2566 1 2 + 404211.06 it 3 
+ 158844096.7 it 5 - 16750593.96 t A + 1912.32 tsi 

- 144 x 4 + 5556.437 t 2 x Sl + 7649.28 tx 2 Sl + 144 Si 2 

+ 418.4064 itxs 1 + 5099.52 itx 4 - 7649.28 itx 2 + 64 x 6 

- 448662.55 t 3 Sl - 101506.4548 it 2 s x + 576 is x x 2 



A 2 =9 + 351839351.0 1 6 - 49079.6020 tV + 33987.4508 tV 
- 5736.96 tSi - 448662.552 t 3 Sl + 5998449.36 tV 
+ 4344308.56 t 5 x + 49254.311 x 3 t 3 + 1605.7344x 3 t 
+ 139.4688 xH - 848937. 199 xt 3 - 1073.5488 tx 
+ 108 x 2 + 78649.7667 1 2 + 48 x 4 + 13134403.21 1 4 
+ 5556.437 t 2 xsi + 7649.28 tx 2 s x + 64 x 6 + 144 Sl 2 

On the other hand, we also give the formulas for q 2 under the condition a — 5, j3 — 1, Si 
0,a = 3.88,c= 1. 



i(3.88 x-30.140928 1) ^3 

q2 ~ e A, 
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+ 20398.08 is tx 3 + 20398.08 is Hx + 22225.748 is 2 t 2 x 
+ 22225.748 it Vs + 3605445.03 it 3 xs - 180 s 2 - 144 s 4 

- 504792.8896 t 2 s x + 11996906.61 i 4 s x + 135949.802 tVs 
+ 203924.713 t 2 x 2 s 2 + 135949.802 s 3 £ 2 :r + 697.3440 s Hx 

+ 1394.6880 sqHx 2 + 697.3440 x 4 s t + 1394.6880 x 3 s 2 t 
+ 3980.3904 x 2 s t + 3980.3904 xs 2 t + 147763. 18 £ 3 s 2 ^ 
+ 147763.15 t 3 s x 2 + 7408.583 is 3 t 2 + 1307537.37 is t 4 

- 127944.917is t 2 + 1802722.57is ^ 3 - 360 s x 
+ 64 s 6 + 5099.52 is H - 7649.28 is 2 t - 576 x 3 s 

A 4 =9 + 64 x 6 + 48 a; 4 + 108 x 2 + 351839241.0 1 6 

+ 13134403.54 1 4 + 78649.7667 1 2 + 5998445.94 tV 
+ 33987.4508 tV - 49079.6020 tV - 848937. 199 t 3 x 
+ 49254.351 tV - 1073.5488 tx + 139.4688 tx 5 
+ 4344312.56 t 5 x + 1605.7344 tx 3 + 49254.311 s 3 t 3 

- 848937.199 s t 3 + 1605.7344 s 3 t + 139.4688 s H 

+ 5998445.94 t 4 s 2 + 33987.4508 s H 2 - 1073.5488 s t 

- 49079.6018 s 2 t 2 + 192 xs 3 + 384 x 5 s + 960 x 4 s 2 

+ 1280 x 3 s 3 + 288 x 2 s 2 + 384 s 5 x + 960 s 4 x 2 + 192 x 3 s 

9 



+ 216 s x + 4344308.56 t 5 s + 108 s + 48 s 4 + 64 s e 



2^3, 



98159.2036 rs x+ 11996906.81 fs x + 135949.802 t z x 6 s 



1.2 2 



3+2, 



+ 203924.713 *Vs + 135949.802 s H z x + 697.3440 s *tx 
+ 1394.6880 s Hx 2 + 697.3440 x 4 s t + 1394.6880 x 3 s H 
+ 4817.2032 x 2 s t + 4817.2032 xs H + 147763. 15 t 3 s 2 x 
+ 147763. 15 t 3 s x 2 
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Figure 1. Higher order rogue waves with zero shifts (fundamental pattern), i.e., (a) n = 
1, a = 10, = 1, a = 0, c = 1, so = 0. (b) n = 2, a = 5, = 1, a = 3.88, c = 1, s = 0, ax = 0. 
(c) n = 3, a = l,/3 = gL, a = 0, c = 1, so = 0, si = 0, S2 = 0. (d) n = 6, a = l,/3 = §g,a = 
0,c = l,so = 0, si = 0,^2 = 0, S3 = 0, S4 = 0,S5 = 0. Here in the figure of \qe\ 2 , to show 
the lower peaks clearly, the whole amplitude of the maximum (central peak) is cut. The 
maximum is 169 (that is (2n + l) 2 )- 
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Figure 2. Single ring structure of \q n \ 2 (left) and the corresponding density plot (right) 



under the condition (a) n = 3, a = 1, /3 



50' 



20 



, c = 1, s = 0, si = 0, s 2 = 1000. The 



trajectory is a pentagon; (b) n = 4, a == 1, j3 = gg, a = 0, c = 1, sq = 0, si = 0, s 2 = 0, S3 = 
1000; (c) n = 5, a = l,/3 = i,a = 0,c = l,s = 0, s x = 0,s 2 = 0, s 3 = 0,s 4 = 4500; (d) 
n = 6,a = l,/3 = ^, a = 0, c = 1, so = 0, si = 0, S2 = 0, S3 = 0, S4 = 0, S5 = 10000. Here 
in the figure of \qe\ 2 , to show the lower peaks clearly, the whole amplitude of the maximum 
(central peak) is cut. The maximum of the inner order 4 rogue wave is 81. 
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Figure 3. Triangular structure of \q n \ 2 (a) n = 2, a 
0,si :- 100; (b) n ^ 3, a - l,/3 - i,< 
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= 5,/3 
l,so 



1, a 
0,si 



3.88, c = 
100, s 2 



n = 4, a = 1, /3 = gg, a = — ^, c = 1, s = 0, si = 100, s 2 = 0, s 3 = 0; (d) n = 6, a 
i, a = 0, c = 1, s = 0, si = 50, s 2 = 0, s 3 = 0, s 4 = 0, s 5 = 0. 



Mo = 
-- 0; (c) 

= l,/3 = 
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Figure 4. The decomposition of |<? n | 2 into double-ring structure (a) n = 4, a = 1,(3 
J_„_n^— 1 e„ — no,— n e„ — 1 nnn B „ — n- fu\ „ — c: „, _ i fl_i„_n^— i Q „ 



^,a = 0,c 



lie ueuuiiipu&iiiun ui |y n | niiu uuuuie-iing suuci 

in 1, s = 0, si = 0, s 2 = 1000, s 3 = 0; (b) n = 5, a = 1, /3 = ^, a = 0, c = 1, s 

0, si = 0, S2 = 0,S3 = 10000, S4 = 0; (c) n = 6, a = 1, /? = ^j, a = 0, c = 1, so = 0, si 
0, s 2 = 0, s 3 = 0, s 4 = 6000, s 5 = 0. 
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Figure 5. Triple-ring structure of \q n \ 2 (a) n = 5, a = l,/3 = X, a = 0,c = 1, sq = 0, s 
0, s 2 = 8000, s 3 = 0, s 4 = 0; (b) n = 6, a = 1, /? 
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0, c = 1, s = 0, si = 0, s 2 = 0, s 3 



1000000, s 4 = 0, s 5 = 0. 
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Figure 6. Ring plus triangular decomposition of \q n \ 2 (a) n = 4, a = l,/3 = 
l, So = o,si = 10, s 2 = 0,s 3 = 10000; (b) n = 5, a = l,/3 = ^,a = 0, c = 1 
50, s 2 = 0, s 3 = 0, s 4 = 10000000; (c) n = 6, a = 1, /3 = i, a = 0, c = 1, s = 0, 
0,s 3 = 0,s 4 = 0,s 5 = 10000000. 
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Figure 7. Fundamental pattern plus ring structure decomposition (a) n = 5, a = l,/3 
^ )G = 0,c = l,s = 0,si = 0,s 2 = 0,s 3 = 80000, s 4 = 20000000; (b) n = 6, a = 1, (3 
^, a = 0, c = 1, so = 0, si = 0, s 2 = 0, s 3 = 0, s 4 = 100000, s 5 = 10000000. 
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Figure 8. Complete decomposition of \qe\ 2 into rings of single peaks a = 1,(3 
0, c = 1, s = 0, si = 0, s 2 = 1000, s 3 = 0, s 4 = 0, s 5 = 100000000. 
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